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ABSTRACT 

We use the results from a constrained, cosmological MHD simulation of the Local Universe to 
predict the radio halo and the gamma-ray flux from the Coma cluster and compare it to current 
observations. The simulated magnetic field within the Coma cluster is the result of turbulent 
amplification of the magnetic field during build-up of the cluster. The magnetic seed field 
originates from star-burst driven, galactic outflows. The synchrotron emission is calculated 
assuming a hadronic model. We follow four approaches with different distributions for the 
cosmic-ray proton (CRp) population within galaxy clusters. The radial profile the radio halo 
can only be reproduced with a radially increasing energy fraction within the cosmic ray proton 
population, reaching >100% of the thermal energy content at w IMpc, e.g. the edge of the 
radio emitting region. Additionally the spectral steepening of the observed radio halo in Coma 
cannot be reproduced, even when accounting for the negative flux from the thermal SZ effect 
at high frequencies. Therefore the hadronic models are disfavored from present analysis. The 
emission of 7-rays expected from our simulated coma is still below the current observational 
limits (by a factor of ~6) but would be detectable in the near future. 
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1 INTRODUCTION 

Galaxy clusters are the largest gravitationally bound objects in the 
Universe. The thermal gas, which forms the dominant component 
in the Intra-Cluster-Medium (ICM), is mixed with magnetic fields 
and relativistic particles, as seen by radio observations that detected 
Mpc- sized diffuse radio sources in a fraction of X-ray luminous 
galaxy clusters in the merging phase (e.g. Feretti 2003 ; Ferrari et al. 
2008). A fraction of the energy dissipated during cluster mergers 
may be channelled into the amplification of the magnetic fields (e.g. 
Dolag et al. 2002 ; Subramanian et al. 2006) and into the accelera- 
tion of relativistic, primary electrons (CRe) and protons (CRp) via 
shocks and turbulence ( |Ensslin et aL]|1998| |Blasi||200l| |Brunetti| 
|& La zarian 2007). CRp have long life-times and remain confined 
within clusters for a Hubble time (e.g. Blasi et al. 2007 and ref. 
therein). Consequently they are expected to be the dominant non- 
thermal particle component in the ICM and to generate secondary 
particles through collisions with thermal protons. Primary and sec- 
ondary particles in the ICM are expected to produce a complex 
emission spectrum from radio to 7-rays (see |Petrosian et al.| 2008 
Brunetti 2008 C assano 2009 for recent reviews). 

Only upper limits to the 7-ray emission from clusters have 
been obtained so far ( Reimer et al. 2003 ; Aharonian et al. 2009), 
although the FERMI telescope will shortly provide the chance to 
obtain the first 7-ray detections of clusters and to put stringent con- 



straints to the energy density of the CRp. Future deep observations 
with high energy Cherenkov arrays are expected to provide com- 
plementary constraints. Most importantly, in a few years the Low 
Frequency Array (LOFAR) and the Long Wavelength Array (LWA) 
will observe clusters at low radio frequencies with the potential to 
discover the bulk of the cluster- scale synchrotron emission in the 
Universe ( [EnBlin & Rottgering|2002l|Brunetti et al.|2008y 

The theoretical picture for 7-ray emission is very complex and 
modern numerical simulations provide an efficient way to obtain 
detailed models of non-thermal emission from clusters to compare 
with present and future observations. Advances in this respect have 
been recently obtained by including aspects of cosmic-ray physics 
into cosmological Lagrangian simulations (Pfrom mer et al.| 2007 
2008), mostly focussing on the acceleration of CRe and CRp at 
shocks and on the production of secondary electrons from such a 
CRp population. In this Letter we investigate the non-thermal emis- 
sion from secondary particles in a Coma-like cluster extracted from 
a cosmological simulation and, for the first time, compare numeri- 
cal predictions and observations. 



2 THE SIMULATION 

We use results from one of the constrained, cosmological MHD 
simulations presented in Don nert et al.| (2009) from which we se- 
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Figure 1. The left panel shows the comparison of the radial profile of the magnetic field in the simulated Coma cluster (black), the scaled version for Model 
3 (blue) and the class models inferred from the observations. The middle panel shows the Energy density fraction of the CRp as function of radius for the 
different models as indicated in the plot. The right pannel shows the radial profile for the radio emission resulting of the different models compared with the 
observed profile. See text for more details. 



lect the simulated counterpart of the Coma cluster. The initial con- 
ditions for a constrained realization of the local Universe were the 
same as used in Mathis et al. (2002). Briefly, the initial conditions 
were obtained based on the the IRAS 1.2-Jy galaxy survey (see 
Dolag et al. 2005 for more details). Its density field was smoothed 
on a scale of 7 Mpc, evolved back in time to z = 50 using the Zel- 
dovich approximation and used as an Gaussian constraint ( Hoffman 
& Ribak 1991 ) for an otherwise random realization of a ACDM 
cosmology (0 M = 0.3, A = 0.7, h = 0.7). The IRAS obser- 
vations constrain a volume of ~ 115 Mpc centered on the Milky 
Way. In the evolved density field, many locally observed galaxy 
clusters can be identified by position and mass. The Coma clus- 
ter in particular can be clearly identified by its global properties 
and gives an excellent match to the observed SZ decrement, espe- 
cially in the non-radiative simulations used in this Letter |Dolag 
|et al.|[2005] >. The similarity in morphology of the simulated and 
observed Coma cluster is coincidental, as this structure is far be- 
low the constrains originally imposed by the IRAS galaxy distri- 
bution. The original initial conditions were extended to include gas 
by splitting dark matter particles into gas and dark matter, obtaining 
particle of masses 6.9 x 10 8 M and 4.4 x 10 9 M respectively. 
The gravitational softening length was set to 10 kpc. 

Our MHD simulation follow the magnetic field through the 
turbulent amplification driven by the structure formation process. 
For the magnetic seed fields a semi-analytic model for galactic 
winds was used. Here we used the result of the 0.1 Dipole sim- 
ulation (Donnert et al. 2009), which gives a reasonable match to 
the observed magnetic field in the Coma cluster. The left panel 
of figure [T] compares the magnetic field profile predicted for the 
Coma cluster from our simulations with models that best repro- 
duce the Rotation Measure observed in five extended sources 
within the Coma cluster (preliminary results from Bonafede et al. 
2009 Bonafede et al. in prep) with a magnetic field radial profile : 

(B(r)) — Bo \ nS nQ V ' > ] an d V i n me range [0.5 -1.5]. Our simu- 
lations are in good agreement with a more centrally peaked profile 
(e.g. r\ « 1.0). This profile gives an average magnetic field over the 
central Mpc 3 of ~ 1.9/iG, which is consistent with the equiparti- 
tion estimate derived from the radio halo emission ( Thierbach et al. 
12003) . 



3 SECONDARY ELECTRONS IN GALAXY CLUSTERS 

The origin of Mpc-scale radio emission in galaxy clusters is still not 
fully understood. Extended and fairly regular diffuse synchrotron 



emission may be produced by secondary electrons injected during 
proton-proton collisions, since the parent relativistic protons can 
diffuse on large scales (e.g. hadronic or seconday models; Den- 
|nison| ( [T980] ); |Blasi & Colafrancesc o (1999)). Alternative models 
assumerelativistic electrons re- accelerated in-situ by MHD turbu- 
lence generated in the ICM during cluster-cluster mergers (e.g. re- 
acceleration models; Brunetti et al. (2001); Petrosian (2001)). 

In this letter we focus on hadronic models that can be im- 
plemented in present simulations. The spectrum of CRp was taken 
to be a power law, N(E) = K p E~ OLp , and the spectrum of the 
secondary electrons resulting from p-p collisions was calculated 
under stationary conditions considering synchrotron and inverse 
Com pton losses ([Brunetti & Blasi|2 005 and ref. therein). Follow- 
ing |Dolag &TS nfilin (2000} we adopted the slope of the spectral 



hadronic models (this also accounts for the logarithmic increase of 
the p-p cross-section with proton energy). Such CRp collisions also 
produce 7-rays from the decay of secondary neutral-pions, which 
are a direct measure of the CRp and of the unavoidable secondary 
electron injection process into the ICM. We calculate the 7-ray flux 
from our simulated Coma cluster by adopting a p — 2.6 and the 
formalism described in Pfrommer & EnBlin (2004). 

The spatial distribution of CRp in our MHD simulations is a 
free parameter, so we follow four approaches chosen to encompass 
a range suggested by theoretical and observational findings, but re- 
stricting to large magnetic fields in the ICM that provide the most 
favourable case for hadronic models. In Model 1 we assume that 
the energy density of the CR protons is a constant fraction of the 
thermal energy density of the ICM. In Model 2 we adopt a radius- 
dependent ratio between CRp and thermal energy density follow- 
ing results from cosmological simulations of CRp acceleration at 
structure formation shocks by Pfrommer et al. ( 2007). In both these 
models, we use the magnetic field strength and spatial distribution 
from our MHD simulations of the Coma cluster. In Model 3 we 
adopt the radius-dependent energy density of CRp as in Model 2 
and an artificial magnetic field by assuming a radial scaling of the 
field strength in the form B oc yfp within the radio emitting region 
(Figure 1). In the last model, Model 4, we used the magnetic field 
from our MHD simulations, but leave the radial profile of the CRp 
energy density in the simulated cluster free to vary to match the 
synchrotron brightness profile measured for the Coma radio halo; 
in this case we limit the CRp energy density so that is stays smaller 
than the thermal energy density of the ICM. 

For all the models the overall normalization of the CRp energy 
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Figure 2. Synthetic radio maps of the simulated Coma cluster (4 Mpc a side). Upper panels are synchrotron maps at 1.4 GHz, lower panels are 7-ray maps 
at E 1 > 100 GeV. The four models displayed from left to right are: constant CRp energy fraction (Model 1), scaled CRp energy fraction (Model 2), scaled 
magnetic field (Model 3) and fitted CRp energy fraction (Model 4). The overall normalization of the CRp energy fraction is chosen to match the observed total 
radio luminosity of Coma. 



density is choosen to match the observed radio luminosity of th e 
Coma radio halo at 1.4 GHz, 7.76 x 10 23 W/Hz |Deiss et al.|l997l. 



4 THE RADIO HALO OF COMA 

The Coma cluster hosts the prototype of giant radio halos (eg. Gio- 
vannini et al. 1993), and in this section we compare the radio prop- 
erties from the simulated Coma cluster with the observed ones. We 
calculate simulated radio images by using a map-making algorithm 
(Dolag et al. 2005), that allows us to project the predicted emission 
of every SPH particle along the line of sight, considering an inte- 
gration depth of ±4Mpc around the center of the simulated Coma 
cluster. The resulting images of radio emission for our models are 
shown in the upper panels of Figure [2] 



4.1 Radial profile and cosmic ray energy budget 

The radial distributions of the fraction between the energy density 
of CRp and that of the thermal gas for Models l^t are reported 
in the central panel of Figure [T] Provided that the magnetic field 
strength in the cluster-central regions is sufficiently large (i.e. > 
5 /iG, as in our cases), Models 1-3 generate enough synchrotron 
luminosity at 1.4 GHz to match that of the Coma halo with rea- 
sonable requests in terms of energy density in the CRp (about 1 to 
10 % of that of the thermal gas). However, a drawback of hadronic 
models discussed in the literature is that the radial distribution of 
the synchrotron emission generated from secondary electrons is ex- 
pected to be much steeper than that observed in radio halos, with 



most of the radio luminosity generated in the cluster-core region 
(e.g. |Brunetti|2004| ). 

The right panel of Figure [T] shows a comparison of the sim- 
ulated and observed radial profile of the radio emission from the 
Coma cluster. In line with previous si mulations ([Dolag & E nBlin 
2000 ) and analytical expectations (e.g. Brunetti 2004]) we find that 
assuming CRp contain a constant fraction of the energy density of 
the thermal pool (Model 7) a profile is produced that is far too steep. 
Flatter radial profiles of the radio emission, and thus fairly extended 
radio halos, are generated only by Models 2 and 3, yet the expected 
radio brightness at 0.5-1 Mpc distance from the cluster-center is 
still >10 times fainter than that observed in the Coma halo. This 
picture also arises from the point-to-point scatter plot between ra- 
dio surface brightness and the X-ray brightness from thermal emis- 
sion (left panel in figure[3]); the solid line is the observed correlation 
obtained by Govoni et al. ( 2001| ). 

To reproduce the radial profile as measured at 1.4 GHz by Deiss 
|et al.|(T997] > out to 0.5-1 Mpc distance, in Model 4 we allow the en- 
ergy density of CRp to vary with distance from the cluster center. 
The resulting radial distribution of the amount of energy in CRp 
relative to that of the thermal pool is shown in the central panel of 
Figure [T] and highlights an expected energetic problem of hadronic 
models: to considerably increase the synchrotron emissivity at large 
distance from cluster center, where the generation of secondary par- 
ticles is inefficient due to the small number density of (target) ther- 
mal protons, a secondary model requires that the energy density of 
CRp be extremely large, comparable with (or even larger than) that 
of the thermal pool. Considering the profile measured with high 
sensitivity by Westerbork observations at 330 MHz (green, dashed 
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Mission 


Veritas (diffuse) 


(r < 1.3 Mpc) 


Egret 


Model 


(> O.lTeV) 


(> O.lTeV) 


(> O.lGeV) 


observed 


< 2.0 x icr 12 




< 3.8 x 10- 8 


Modell 


6.3 x 10- 14 


6.6 x icr 14 


1.1 x 10- 9 


Model2 


6.3 x icr 14 


7.8 x icr 14 


1.2 x 10- 9 


ModeU 


6.4 x icr 14 


8.0 x 10- 14 


1.3 x 10- 9 


ModeU 


1.6 x icr 13 


3.7 x 10- 13 


6.0 x 10- 9 



Table 1. Expected limits of 7-ray emission from the simulated Coma clus- 
ter in 7/s/cm 2 , excluding the substructure on the lower right. For the three 
models of CR distribution we show the limits in energy ranges, correspond- 
ing to recent experiments. 



line in the right panel of figure [T] Govoni et al. 2001 ) makes this 
point even more problematic. 



4.2 The Spectrum 

The radio spectrum of the Coma halo shows a steepening at higher 
frequencies (Thierb ach et al.|[2003| ) that has been interpreted as 
a signature of stochastic reacceleration of the emitting electrons 
( |Schlickeiser et al.|1987l|Brunetti et al.|2Q01| . On the other hand, 
it has also been argued that the steepening is due to incorrect esti- 
mation of the flux of the radio halo at the highest frequencies due 
to the Sunyaev-Zeldovich (SZ) decrement (e.g. EnBlin 2002). The 
latter results in a flux reduction of the cosmic microwave back- 
ground photons in the cluster region by Compton scattering with 
the cluster-thermal electrons, which diminishes the radio flux at 
higher frequencies, although other authors conclude that this is 
not sufficient to explain the spectral steepening in the Coma halo 
( |Reimer et aL|2004[[Brunetti|2004| >. 

Our numerical simulations can be used to investigate SZ- 
decrement effect on the synchrotron spectral properties. Figure [3] 
(central panel) shows the radio and inverse Compton spectrum from 
the SZ effect inside the cluster gas. The red line is the prediction 
of the total flux from our models, where the SZ signal is extracted 
from the region of the radio halo at 1.4 GHz, ^500 kpc in radius. 
The deviation from a pure power law at higher frequencies is due 
to the SZ flux decrement (the dotted line marks the corresponding 
power law with v 1 -25 ), while the flattening at lower frequencies is 
due to the energy-dependent cross-section of p-p collisions. We plot 
the absolute of the IC flux, although it corresponds to negative flux 
at frequencies smaller than « 2 x 10 5 MHz. In black we show ob- 
servations from |Thierbach et al.| ( [2003| ) and |Battistelli et al.| ( [2002] ) 
of the Coma cluster in the radio and microwave band, respectively. 
We also evaluated the fluxes integrated over the relevant sizes ac- 
cording to the individual observations (red symbols). The observed 
shape of the radio spectrum, with the steepening at 2 x 10 3 MHz, 
cannot be explained by our secondary models, even including the 
SZ decrement for these frequencies, although the simulated SZ sig- 
nal almost perfectly fits observations in the microwave regime. The 
blue dashed line shows the inverse Compton decrement, assum- 
ing an isothermal beta model for the thermal gas of Coma with 
emission region of 5 Mpc radius (as done in EnBlin (2002)), which 
would be needed to cause the steepening at 5 GHz. This size is 
about one order of magnitude more extended than the radio emis- 
sion. Additionally we show in green the expected spectrum from a 
stochastic electron- acceleration model ( Schlickeiser et al.|1987) . 



5 7-RAY SPECTRUM AND LIMITS 

Maps of the predicted 7-ray emission of our simulated Coma clus- 
ter are shown in the lower panels of Figure [2] In the right panel of 
Figure[3]we show the differential and integral 7-ray flux as function 
of energy of the photons. Here we also included the observed lim- 
its from VERITAS and EGRET ( |Reimer et d.|2003||PerEns|2008) , 
see Table [T] 

We find that the 7-ray emission produced in the case of Model 
4, to allow a reasonable match with the Deiss et al. profile of the 
Coma halo (at least at < 0.7-1 Mpc distance from cluster center), 
is a factor ^6 below present upper limits. As well, matching the 
radial profile of the Coma halo as measured by high sensitivity ob- 
servations at 330 MHz would require an even larger energy budget 
of CRp with respect to that in Model 4 and consequently a larger 7- 
ray emission. Thus as soon as future 7-ray observations will reach 
sensitivities slighly deeper than present upper limits they will allow 
complementary tests to the hadronic origin of the radio emitting 
electrons in the Coma cluster. For instance, the expected sensitivity 
of the Fermi Gamma Ray Telescope will be sufficient to constrain 
Model 4. 



6 CONCLUSIONS 

Based on a constrained, cosmological MHD simulation of the local 
universe we investigated the predicted properties of the radio halo 
of the Coma cluster within the framework of the hadronic model. 
We follow four approaches, chosen to span the reasonable range 
suggested by theoretical and observational findings, and focusing 
on the case of high magnetic field values that represents the most 
favorable way for hadronic models. 
Our main conclusions are: 

• In agreement with previous findings, hadronic models may 
produce the synchrotron radio luminosity of the Coma halo with 
energy densities of CRp between 1-10% of the thermal ICM. How- 
ever the radial brightness profile of the generated synchrotron emis- 
sion is much steeper than seen in observations of the Coma halo, 
consequently the simulated radio halos come out much smaller than 
the observed one. This also leads to a slope in the thermal X-ray 
vs. radio brightness point-to-point correlation that is significantly 
steeper than the observed one. 

• The observed flat radial brightness profile and the fairly large 
extent of the observed radio halo can only be obtained with 
hadronic models by strongly increasing the CRp energy density 
outside the core of the cluster. In this case, however, the resulting 
CRp energy density at w IMpc (e.g. the rim of the observed radio 
halo) would equal (or exceed) that of the thermal ICM. 

• Our simulated Coma cluster matches almost perfectly the ob- 
served SZ decrement. But this SZ decrement is still not enough 
to explain the spectral steepening observed in the spectrum of the 
radio halo at large frequencies; contrary to previous claims in liter- 
ature. 

In summary, we find that the hadronic model is disfavored by cur- 
rent observations of the Coma cluster. In principle the energy prob- 
lem might be alleviated by if the magnetic field strength were al- 
most constant with cluster radius up to the rim of the observed radio 
halo (see Pfrommer & EnBlin 2004), however such models are not 
among the best-fitting ones inferred from present RM observations. 
At the same time we have shown that hadronic models cannot ex- 
plain the spectral steepening of the Coma radio halo: regardless of 
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Figure 3. The left panel shows the Radio Flux of the simulated Coma cluster at 1.4 GHz versus X-ray surface brightness. We plot the constant model in 
black, the model with scaled CRp in red, the one with scaled MF and CRp in blue and the one with fitted CRp in magenta. The solid line is taken from the 
observations by Govoni et al. ( 2001). The middle panel shows the observed spectrum of the Coma cluster (synchrotron emission: black diamonds, Thierbach 
et al. ( 2003); inverse Compton: black triangles, Battistelli et al. (2002)). Red lines are from the simulation and the red symbols are evaluated within the same 
aperture as the according observations. See text for more details. The right panel shows the expected 7-ray spectrum from the simulated Coma cluster for 
the our different models for the CRp distribution and the integrated 7-ray spectrum (within r < 1.3 Mpc) with the current upper limits from EGRET and 
VERITAS (Reimer et al.T200^|Perkins|2008> . 



the model assumptions, an SZ effect matching observations pro- 
duces a negligible SZ decrement at 2.7 and 5 GHz in the region of 
the halo. 

The 7-ray flux generated by the hadronic interaction in case 
of Model 4, which alone matches the observed radial profile is only 
a factor ^6 below the current limits. Considering that in case of 
lower magnetic fields (and/or assuming the more extended profile 
of the Coma halo from high sensitivity 330 MHz observations) the 
expected 7-ray flux should significantly increase, we conclude that 
incoming 7-ray observations will shortly provide additional con- 
straints to the models. 
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